Astronomy & Astrophysics manuscript no. paper 


© ESO 2009 


July 14, 2009 





Measuring cosmic shear with the ring statistics 



T. Eifler 1 , P. Schneider 1 and E. Krause 2 ' 1 



o 
o 



1) Argelander-Institut fur Astronomie, Universitat Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 

2) California Insitute of Technology, M/C 350-17, Pasadena, California 91125, USA 



ABSTRACT 



o 

u 

43 

Or 
i 

o 



> 
o 

m 

(N 

o 

O 



Context. Commonly used methods to decompose E- and B-modes in cosmic shear, namely the aperture mass dispersion and the E/B-mode 
shear correlation function, suffer from incomplete knowledge of the two-point correlation function (2PCF) on very small and/or very large 
scales. The ring statistics, the most recently developed cosmic shear measure, improves on this issue and is able to decompose E- and B-modes 
using a 2PCF measured on a finite interval. 

Aims. First, we improve on the ring statistics' filter function with respect to the signal-to-noise ratio. Second, we examine the ability of the 
ring statistics to constrain cosmology and compare the results to cosmological constraints obtained with the aperture mass dispersion. Third, 
we use the ring statistics to measure a cosmic shear signal from CFHTLS (Canada-France-Hawaii Telescope Legacy Survey) data. 
Methods. We consider a scale-dependent filter function for the ring statistics which improves its signal-to-noise ratio. To examine the 
information content of the ring statistics we employ ray-tracing simulations and develop an expression of the ring statistics' covariance in 
terms of a 2PCF covariance. We perform a likelihood analysis with simulated data for the ring statistics in the Q m -<r 8 parameter space and 
compare the information content of ring statistics and aperture mass dispersion. Regarding our third aim, we use the 2PCF of the latest 
CFHTLS analysis to calculate the ring statistics and its error bars. 

Results. Although the scale-dependent filter function improves the S/N ratio of the ring statistics, the S/N ratio of the aperture mass dispersion 
is higher. In addition, we show that there exist filter functions which decompose E- and B-modes using a finite range of 2PCFs (£B-statistics) 
and have higher S/N ratio than the ring statistics. However, we find that data points of the latter are significantly less correlated than 
data points of the aperture mass dispersion and the £6-statistics. As a consequence the ring statistics is an ideal tool to identify remain- 
ing systematics accurately as a function of angular scale. We use the ring statistics to measure a E- and B-mode shear signal from CFHTLS data. 

Key words, cosmology: theory - gravitational lensing - large-scale structure of the Universe - methods: statistical 



1. Introduction 



Cosmic shear was first detected in 2000 dBacon et alj 2000t 



Kaiser et al. 2000t van Waerbeke et al. 2000I ; IWittman et al 



20001) and has progressed to a valu able source of cosmologi 



arate the cosmic shear signal from intrinsic alignment contam- 
inations if redshift information is available. The strength of B- 
modes coming from these effects are examined t hrough numer- 



ical s i mulations; al though th e results di ffer (e.g. lHeavens et al 



cal information. Latest results (e.g.. | va n Waerbeke et alll2005t 
Semboloni et alj 2006 : Hoekstra et alJbOOd ISchrabback et al 



2007 : iHetterscheidt et al.ll2007t iMassev et al.ll2007fj; iFu et al 
20081) already indicate its great potential to constrain cosmo- 
logical parameters, which will be enhanced by large upcoming 
surveys like Pan-STARRS, KIDS, DES or Euclid. 
An important step in a cosmic shear analysis is the decom- 
position into E- and B-modes, where, to leading order, grav- 
itational lensing only creates E-modes. In principle, B-modes 
can arise from the limited validity o f the Born approximation 



l200d;ICrittenden et alj|200ll;|jindl2002l) . the observed B-mode 
amplitude is higher than expect ed from the foregoing expla- 
nations. Shape-shear correlation ( IHirata & Seliakll2004l) is an- 
other astrophysical contamination w hich can cause B-modes. 
Joachimi & Schneider! (|2008, 120091) show how to exclude the 



dJain et alj|200d; iHilbert et al 



2008) or redshift source clus- 



tering dSchneider et alj l2002bl) . Another possible source are 



astrophysical co ntaminations such as int rinsic alignment of 
source galaxies; iKing & Schneider! d2003l) show how to sep- 



contaminated scales, again using redshift information. 
Most likely, B-modes indicate remaining systematics in the ob- 
servations and data analysis, in particular they can result from 
an insufficient PSF-correction. The Shear TEsting Program 
(STEP) has significantly improved on this issue (for latest re- 
sults see Heymans et al. 20061 : Massey et al. 2007al) : still the 
accuracy of the ellipticity measurements must be improved fur- 
ther to meet the requirements of precision cosmology. 
The identification of remaining systematics (B-modes) will be 
important especially for future surveys, where the statistical er- 
rors will be significantly smaller. Therefore, decomposing the 
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shear field into E- and B-modes must not be affected from in- 
herent deficits. The most commonly used methods for an E- 
and B-mode decomposition, the aperture mass dispersion and 
the E/B-mode shear correlation function, require the shear two- 
point correlation (2PCF from now on) to be known down to 
arbitrary small or up to arbitrary large angular separations, re- 
spectively. This is not possible in practice; as a consequence the 
corresponding methods do not separate E- and B-modes prop- 
erly on all angular scales. A detai led analysis of this issue can 
be found in lKilbinger et all d2006h (hereafter KSE06). 



Most cosmic s hear analyses, e.g. iMassev et al.l (I2007bb and 
Fu et al l d2008l) (hereafter FSH08), simulate 2PCFs from a the- 
oretical model of P K to account for the scales on which the 
2PCF cannot be obtained from the data. This ansatz is problem- 
atic, since one explicitly assumes that the corresponding scales 
are free of B-modes. In addition, the assumed cosmology in the 
theoretical power spe ctrum can bias the resul ts. 



The ring statistics ([Schneider & Kilbingerl 120071 hereafter 
SK07) provides a new method to perform an E-/B-mode de- 
composition using a 2PCF measured over a finite angular range 
[i?min; #max]- In this paper we examine the ring statistics in de- 
tail; more precisely we improve the ring statistics' filter func- 
tion with respect to its S/N ratio and examine its ability to con- 
strain cosmological parameters. Furthermore, we construct a 
filter functions which has higher S/N ratio than the ring statis- 
tics but still decomposes E/B-modes with a 2PCF measured 
over a finite range. We will refer to this as fiB-statistics. 
Due to the fact that the ring statistics' data points show sig- 
nificantly lower correlation than data points of the aperture 
mass dispersion and the fiB-statistics, it provides an ideal tool 
to identify remaining systematics in cosmic shear surveys de- 
pending on the angular scale. We employ the ring statistics to 
identify B-modes in the CFHTLS survey. 
The paper is structured as follows: In Sect. |2]we start with the 
basics of second-order cosmic shear measures, followed by the 
main concepts of the ring statistics in Sect. [3] We derive a for- 
mula to calculate the ring statistics' covariance from a 2PCF 
covariance in Sect. [5] and also compare the correlation coef- 
ficients of ring statistics, aperture mass dispersion, and EB- 
statistics in this section. In the same section we examine the 
S/N ratio of the ring statistics and compare it to the other mea- 
sures. More interesting than the S/N ratio however, is the abil- 
ity of a measure to constrain cosmology. This, in addition to 
the S/N ratio, depends on the correlation of the individual data 
points. In order to quantify this accurately, we perform a like- 
lihood analysis in Sect. [6] for the ring statistics, aperture mass 
dispersion and £B-statistics using data from ray-tracing simu- 
lations. The results of our analysis of CFHTLS data using the 
ring statistics are presented in Sect. [7] followed by our conclu- 
sions in Sect. [8] 

2. Two-point statistics of cosmic shear 

In this section we briefly review the basics of second- 
order cosmic shear measur es. For more details on th i s topi c 



the reader is 



Schneider et a 
Munshi et alJ d2008h" 



r eferred to 
j2002albh : 



Bartelmann & Schneidet] d200ll) ; 
van Waerbeke & Mellieri J2003I) : 



To measure the shear signal we define & as the connecting 
vector of two galaxy centers and specify tangential and cross- 
component of the shear y as 



= -Re( r e- 2i ^) 



and 



y> 



= -Im(ye- 2i ^) , 



(1) 



where <p is the polar angle of The 2PCFs depend only on the 
absolute value of They are defined in terms of the shear and 



can be related to the power spectra Pe and Pb dSchneider et al 
2002bl) 



t±m = <7t7t>(#) ± <rxrx>(#) 

Zn 



(2) 
(3) 



with J„ denoting the «-th order Bessel-function. 
Starting from the 2PCF as the basic observable quantity, there 
exist several methods to decompose E-modes and B-modes, 
such as the E/B-mode shear correlation function or the aperture 



mass d ispersion (e.g. ICrittenden et al.l l2002t ISchneider et al 
2002b). The latter can be calculated as 



1 f 2H d#i? 

The filter functions read 
6(2 - 15x 2 ) 



f + (0)r + |£)±f_(0)r-(£ 



(4) 



T + (x) = 



T-(x) 



1 arcsin 

n 



(I) 



IOOtt 



(120 



+2320x 2 - 754/ + 132x 6 - 9x 8 )} H(2 - jc) , 



192 
35^' 



je 3 1 1 



2\7/2 



H(2 - x) ., 



(5) 



(6) 



with H being the Heaviside step function. Decomposing E- and 
B-modes with the either the aperture mass dispersion or the 
E/B-mode shear correlation function requires that the 2PCF is 
either measured down to arbitrary small or large angular sep- 
aration, respectively. For further details on this problem the 
reader is referred to KSE06. 

3. The ring statistics 

To circumvent the aforementioned difficulties SK07 in- 
troduced the ring statistics whose second-order moments 
({'R'Re) , (7?9?b)) decompose E- and B-modes properly using 
2PCFs measured on a finite interval [# m i n ; #max]- The quan- 
tity CRKe) can be interpreted as the correlator of the shear 
measured from galaxy pairs which are located inside two con- 
centric rings (see Fig. [T]). Their annuli are chosen as follows: 
<fi < Q\ < £> for the first ring and £3 < 62 < (4 for the sec- 
ond. The rings are non-overlapping, i.e. £j < gj if i < j. The 
argument of the rings statistics is named *P = (2 + Ca and only 
2PCFs with & < m enter in the calculation of {"R'Re}. In addi- 
tion, the ring statistics depends on a parameter r\ quantifying 
the separation between outer and inner ring, i.e. rj/W = (3- 
In order to calculate the ring statistics properly from a set of 
2PCFs within [# m j n ; #max] it is required that *P does not exceed 
i? max and that i? m i n / > T < 77 < 1. 

Following the derivation of SK07 the E- and B-mode decom- 



T. Eifler, P. Schneider and E. Krause: Measuring cosmic shear with the ring statistics 



3 




Fig. 1. This figure illustrates the basic idea of the ring statis- 
tics and how it can be obtained from the 2PCF of cosmic shear. 
We measure the 2PCF of each galaxy in the inner ring with 
all galaxies in the outer ring. For a given argument of the ring 
statistics *P, the angular separation of the required 2PCFs ex- 
tends over rfY <§ < X V. The meaning of 77 and its possible val- 
ues are further explained in the text. The ring statistics is then 
calculated as an integral over the 2PCF with the filter functions 

position of the ring statistics can be obtained from the 2PCF 

as 



<^B>CP) 



r' r d# 

[?+(#) Zf(0,J7) +£-(#) Z~(0,J7)] , (7) 
Zf(#, n) - Z-(0, n)] ■ (8) 

J„<v 2# 



The functions Z ± are defined in SK07; we plot them in Fig. [2] 
for four different 77, i.e. ^n/T = 0.00151,0.1,0.4,0.7. 
Similar to the case of the aperture mass dispersion, {'R'Re) can 
be related to the E-mode power spectrum. Inserting Eq. (0, 
into Eq. (0 gives 



with 



■r 
-r 



~2n 



(9) 



dtf 
2# 



[j o (€ff)z + (0,Ti) + Mmz-(&,n)] ■ (10) 



When calculating {'R'Re) for different arguments W, we distin- 
guish two cases for 77. It can be fixed to a specific value or it 
can vary according to *F, in particular 77 = !? m in/ v f. We will 
refer to the latter case as a scale-dependent 77. Here, the lower 
limit in the integrals of Eqs. (0 and ([8]l is equal to m j n which 
implies that all 2PCFs in the interval [0 m i n ; W] are included in 



the calculation. The choice of 77 = 7? m i n / v I' should give a higher 
S/N ratio compared to a fixed 77 for the reason that more galaxy 
pairs are included which reduces the statistical noise. In SK07 
the authors hold 77 fixed; in order to obtain a high signal this 
implies that 77 must be chosen as small as possible. 
Choosing a fixed 77 has a second disadvantage. The lower limit 
in the integrals Eqs. © and <[8J cannot be smaller than m i n , i.e. 
77 > T^min- Vice versa, this implies that > tfmin/'?- Fixing 77 
to a small value (in order to increase the S/N ratio) implies that 
\P is restricted to larger scales. This trade-off between S/N ratio 
and small-scale sensitivity can be overcome when relaxing the 
condition of a fixed 77. 

4. General E/B-mode decomposition on a finite 
interval 

The ring statistics described in the last section is the special 
case of a general E/B-mode decomposition. According to SK07 
this general /iB-statistics can be defined as 



1 

2 j„ 
1 



d00[£ + (0)7 + (0) + £.(0)71(0)] , 



B=-\ d00[f + (0)7 + (0)- £_(0)7_(0)] 
^ Jo 



(11) 
(12) 



To provide a clean separation of E- and B-modes using a 2PCF 
measured over a finite interval, the following conditions must 
be fulfilled (see SK07 for the exact derivation). Starting from 
an arbitrary function T + {ff), which is zero outside the interval 
[#min; #maxL the constraints 



f 



d007 + (0) = = 



d0# 3 71(#) 



(13) 



must hold. For a so constructed filter function T + (§) a corre- 
sponding filter function 7_(0) can be calculated as 



7_(0) = 7 + (0) + 4 £ dfl|r + (fl)[l-3Q 



(14) 



Conversely, one can construct T + for a given 7_. 

The expressions for T + and T used in this paper are given in 

the Appendix. We calculate the /iB-statistics according to Eq. 

(fTTT i and compare the results to the ring statistics. Note that this 

fiB-statistics can be optimized, e.g., with respect to its S/N ratio 

or its ability to constrain cos mology. For more deta ils on this 

topic the reader is referred to lFu & Kilbingerld2009l) . 

In this paper, the EB-statistics is calculated as a function of *P. 

Similar to the ring statistics, *P denotes the maximum angular 

scale of 2PCFs which enter in the calculation of £(*?). 

5. Covariance and signal-to-noise ratio 

For our further analysis we have to derive a formula to calcu- 
late the covariance of ring statistics and ffB-st atistics. A corre- 



spondi ng expression for (M^) reads (see e.g. ISchneider et al 



2002b). 



1 ^, <^ A0/A0; 

c M (o k ,e,)) = j2-j£j 

'=1 j=i 
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Fig. 2. This plot shows the filter functions Z+ (left ) and Z (right) depending on # for four different choices of rj: # m in/ v P = 
0.00151 (solid), 0.1 (dashed), 0.4 (dotted), 0.7 (dotted dashed). 



2 Mt) r-fD^GW 



(15) 



with C m „(i?,, denoting the 2PCF covariance. Here, the upper 
limits / and J are chosen such that < 2ft and §j < 28/. The 
ring statistics' covariance is defined as 

CrCF*. Y£> = - <^e>CP0 <^ E > W) , (16) 

where ^ 2 denotes the estimator of the ring statistics. To cal- 
culate this estimator from a binned 2PCF data vector with bin 
width A#, we replace the integrals in Eq. (0 by a sum over the 
bins 

SePP) = \ Yj ^ Z +^<' ^ + Z -^' ^] ' < 17) 

!=1 ' 

with f±(#j) denoting the estimator of the ;-th 2PCF bin. The 
upper limit / in Eq. ( TP71 ) denotes the bin up to which < *P. 
Inserting Eq. ( TPTI i into Eq. ( fTBT ) we derive 



(18) 



2 ZJftu W k ) Z n (#j, %) CM, Oj) 

m,n=+- 

where / and J denote the bins up to which &, < *Pt (§j < ¥/) 
holds. 

Similarly a covariance for the general £B-statistics can be cal- 
culated as 

/ j 

i=l ;=1 



2 r m (#;, e,) r„ ft) c™^,-, 



(19) 



5.7. Correlation matrices 

In order to illustrate the correlation between the individual data 
points we calculate the correlation matrix R for {'RfR^), E, and 
<M 2 p ) from the corresponding covariance matrix. For C being 
the covariance of either (KRe), E, or <M 2 p ) the correlation co- 
efficients are defined as 



Ry = 



(20) 



The covariances are calculated from a 2PCF ray-tracing co- 
variance via Eqs. ( fT5l ), ( fT8l , and ( fT9l ), respectively; finally the 
correlation matrix is obtained via Eq. ( |20l . The ray-tracing sim- 
ulations (175 realizations) have the following underlying cos- 
mology: Q m = 0.27, C2 A = 0.73, cr 8 = 0.78, h = 0.73, Q b = 
0.044, n s = 1 .0. From now on we refer to this cosmological pa- 
rameter set as our fiducial cosmological model jify. Survey pa- 
rameters which enter in the calculation read as follows: galaxy 
density n ga i = 25/arcmin 2 , survey area A — 36 deg 2 , and in- 
trinsic ellipticity noise cr f = 0.38. The survey parameters differ 
slightly from those of the covariance used in the latest CFHTLS 
survey; FSH08 use A = 34.2 deg 2 , « ga i = 13.3/arcmin 2 , and 
cr e = 0.42. 

The covariance matrices have a different angular range corre- 
sponding to the data vectors of CRRe), E, and (M 2 p ), which we 
define as 



<7?7?e> = [<««E>( , Pi),...,<W?E)( l P„)] t , 

E = [EQFO,...,ECFjf , 
<M 2 p > = [<M 2 p >(ft),..,<M 2 p >(ft0f. 



(21) 
(22) 

(23) 



Whereas <7?7? E > and E extend from 1' < "P < 460', <M 2 p > 
extends from 6' < 6 < 230'. The different maximum angular 
separation of the aperture mass dispersion result from the fact 




Fig. 3. This figure shows the correlation matrices of (KRe) (left), E {middle), and (M 2 p ) (right) derived from ray-tracing 2PCF 
covariance matrix. In each panel the n-th contour line (starting with n — 1 close to the diagonal) marks values of (0.8)". 



that CRUe) (T) and £(¥) contain information on the 2PCF with 
# < "P, whereas {M 2 p )(ff) contains information on the 2PCF 
up to ■& < 28. The lower limit of 6' was chosen to circumvent 
the problem of E/B-mode mixing for the (AfL) covariance. The 
range of the original 2PCF ray-tracing covariance extends from 
0.5 < § < 460'. Below 6' it is not possible to calculate the 
(M 2 p ) covariance properly. 

Figure [3] shows the correlation matrices of the ring statistics 
(left), the SB-statistics (middle), and of the aperture mass dis- 
persion (right). Starting from the diagonal, where /?,-,- = 1, the 
n-th contour line corresponds to values of 0.8". It is clearly no- 
ticeable that data points of the ring statistics are significantly 
less correlated than those of the aperture mass dispersion and 
the SB-statistics. 

The boxy contours in Fig. [3] result from the small number of 
bins we choose in the covariances. The reason for this is that 
the ray-tracing covariance is an estimated quantity; its inverse, 
needed for the likelihood analysis in Sect. [6] is in general af- 
fected from numerical artifacts. These artifacts become more 
severe in case of a high dimension matrix. In order to guar- 
antee a stable inversion process we choose a small number of 
bins. 



5.2. Signal-to-noise ratio 



We now use the above derived covariances to quantify the S/N 
ratio of the ring statistics, SB-statistics and compare both to 
that of the aperture mass dispersion. 

We calculate a set of 2PCFs via Eq. (O for an angular range 
similar to that of the ray-tracing simulations (see Sect l5.lt , i.e. 
i? e [015,460']. The required shear power spectra P E are ob- 
tained from the density power spectra P s employing Limber's 
equation. As underlying cosmology we choose our fiducial 
model (see Sect. [5]). The power spectrum Pg is calculated from 
an initial Harrison-Zeldovich power spectrum (P/j(k ) oc fc* s ) 
with the transfer function from lEfstathiou et"aD ( 1992|). For the 
non-linear evolution we use the fitting formula of ISmith et al 



(|2003l) . In the calculation of Be we ch oose a redshift d i stribu- 
tion of source galaxies similar to that of lBenjamin et al.l (120071) 



n(z) 



P 



z T((l+a)/B)\z 



— I exp 



(24) 



with a = 0.836, B = 3.425, z = 1-171. 

From this set of 2PCFs we calculate data vectors of CR!Re), 
E, and (M 2 ) according to Eqs. (O and ( fTTT ). and (0]), respec- 
tively. The angular range of these data vectors are chosen sim- 
ilar to the range of the corresponding covariances (Sect. 15. It , 
i.e. 0:5 < ¥ < 460:0 for <«« E > and E, and 6'.0 < 6 < 230:0 
for (M 2 ). The S/N ratio is calculated as 



S/N = 



WTO 



[CkOM-)] 



1/2 



and S/N 



{M 2 dp )(6i) 



1/2 



(25) 



The results are illustrated in Fig. |4] We compare the ring 
statistics for with scale-dependent 77 and 77 = 0.1 to the EB- 
statistics and the aperture mass dispersion. The figure shows 
the anticipated behavior (Sect. [3); the ring statistics with scale- 
dependent n gives a larger S/N ratio compared to the case where 
Tj is fixed. In addition, it can be measured down to arbitrary 
small values of *P (above # m j n ), which is not possible when 
choosing a fixed 77. For the case considered here, i.e. # m i n = 0.5, 
the choice of 77 = 0.1 already limits the range of to scales 
> 5'; decreasing 77 further in order to increase the S/N ratio will 
limit CRKe) to larger V P. 

When comparing the ring statistics to the aperture mass dis- 
persion, we find that the ring statistics' signal is lower. Even 
with the scale-dependent filter function the S/N ratio of the 
ring statistics is on average by a factor of « 2 smaller than the 
S/N ratio of the aperture mass dispersion. This difference can 
be explained when comparing the filte r functions of CRRy) an d 
(Ml p ), Z ± (Fig.EJ and T ± (e.g. Fig. 1 in lSchneider et al.l2002bh . 
respectively. The Z-functions have two roots at their boundaries 
whereas the TVfunction becomes particularly large for small x. 
However, we point out that the S/N ratio does not solely deter- 
mine the ability of a measure to constrain cosmology, but one 
has to account for the fact that the data points of the CR!Re) are 
less correlated than those of (M 2 p ). For a full comparison of the 
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Fig. 4. The S/N ratio of the ring 
statistics (for r\ = 01 5/4* and 
r\ = 0.1), the £B-statistics, and the 
aperture mass dispersion calculated 
from a set of theoretical 2PCFs with 
!? e [0:5;460']. The different an- 
gular range of the measures is ex- 
plained in the text. 



information content we examine both measures in a likelihood 
analysis. 

Compared to the ring statistics the S/N ratio of the /sB-statistics 
is significantly larger on all scales, which again can be ex- 
plained by the fact that the filter function of the ZiB-statistics 
does not have roots at their boundaries. Compared to the aper- 
ture mass dispersion, the EB-statistics' S/N ratio is slightly 
lower. However, we point out that the EB-filtei function, we 
chose here, is a simple second-order polynomial. We will 
present an extended analysis of this general EB-filtei functions 
in a future paper. 

6. Comparison of the information content of CRRe) 

and (Ml v ) 

We now perform a likelihood analysis in the Q m vs. erg, pa- 
rameter space in order to compare the ability of CRKe), E, 
and {M%) to constrain cosmological parameters. We calculate 
2PCF data vectors for various combinations of <x 8 e [0.4; 1 .4] 
and Q. m e [0.01; 1.0], therefrom derive the data vectors of 
CRKe), E, and (MjL) and test these against the corresponding 6 
data vectors obtained from our fiducial model (Sect. 15. 11 1. We 
assume that all data vectors are normally distributed in param- 
eter space and calculate the posterior likelihood according to 
Bayes theorem. Our likelihood function p{d\7i) then reads 

exp \-\ Ud{n) - d(7t M )f C- 1 (d(7r) - d(n M )))] 
P(d\n) = j ^ ,(26) 

(27T)"/ 2 |C|5 

where d must be replaced by the considered data vector, either 
<7W? E > (Eq.EB , (Ml p ) (Eq.EUl, or E (Eq.|23]>. 
To illustrate the information content we calculate the so-called 
credible regions, where the true parameter is located with a 



probability of 68%, 95%, 99,9%, respectively. In addition, we 
quantify the size of these credible regions through the deter- 
minant of the second-order moment of the posterior likelihood 
(see e.g. lEifleret ai1l2008al) 



fly = j A 2 np{n\0 (m - afXff, - n]) , 



(27) 



where n, are the varied parameters, n f are the parameter of the 
fiducial model (i = 1,2, corresponding to O m and cr g ). The 
square root of the determinant is given by 



q = y j\Q~\ = ^£n£ 2 2-£i 2 , 



(28) 



and it can be considered as our figure-of-merit quantity. 
Smaller credible regions in parameter space correspond to a 
smaller value of q. In this paper all g's are given in units of 

io- 4 . 

For the likelihood analysis in this section we employ the ray- 
tracing covariances and choose the angular range of the data 
vectors correspondingly (Sect. 15. Il l, i.e. *P e [l';460'] and 
G [6'; 230']- We further assume a flat prior probability with 
cutoffs, which means p(7t) is constant for all parameters inside 
a fixed interval (Q m e [0.01; 1.0], cr 8 e [0.4; 1.4]) and p{n) = 
else. 

As we obtain our covariance from ray-tracing simulations we 
automatically account for the non-Gaussianity of the shear 
field, however we neglect the cosmology depend ence of the co- 
variance (for more details see Eifl er et al . 2008b). Furthermore, 
we account for the bias which occurs during the inversion of the 
ray-trac ing covariance by a pplying the correction factor out- 
lined in lHartlap etail d2007t) . 

The upper row of Fig.|5]shows the result of the likelihood anal- 
ysis for the ring statistics. We consider 3 cases: First, CRKe) 
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Fig. 5. The 68%-, 95%-, 99.9%-contours of 
the likelihood analysis using the ring statis- 
tics, the £B-statistics, and the aperture mass 
dispersion. We compare 5 different cases, 
namely in the upper row: 77 = ¥/#„,(„ for 
fmin = 1' (left), and for = 6' (mid- 
dle), and CRRe) with rj = 0.1 (right). In 
the lower row we see <M^ p ) (left), and the 
ES-statistics for T,,,^ = 1' (right). The data 
vectors are calculated analytically from a 
power spectrum; the covariance is obtained 
from ray-tracing simulations. The filled cir- 
cle marks our fiducial cosmology. 



with T] = tfmin/Y and ¥ € [l';460'] (left). Second, CR!R E ) with 
r, = # min /*F and *F e [6'; 460'] (middle). Third, CR!R E ) with 
77 = 0.1 and e [l';460'] (right). The lower row shows a 
similar analysis for (M^ p ) with an angular range 9 e [6'; 230'] 
(left) and the £B-statistics for *P 6 [l';460'] (right). The black, 
filled circle indicates the fiducial cosmology; the contours cor- 
respond to the aforementioned credible regions. In addition we 
quantify the information content by the values of q, defined in 
Eq. ( l28l i. which are summarized in TableQ] 
The ring statistics with 77 = i? m i n / v f is a clear improvement over 
CR!Re) with 77 = 0.1 which can be explained by the larger S/N 
ratio of the first compared to the second. Considering the ring 
statistics with scale-dependent 77, we find that adding informa- 
tion below 6' increases the information content of CRKe), such 
that it gives tighter constraints than the (M~ p ) data vector. The 
strength of this gain in information can be explained by the 
small correlation of ring statistics' data points. 
In our analysis it was not possible to calculate (0L) for 6 < 6' 
due to the aforementioned E/B-mode mixing, however this can 
change if the 2PCF is measured on smaller angular scales. For 
this case we expect the improvement of ring statistics over the 
aperture mass dispersion to be even more significant. Due to 
the lower correlation of the ring statistics' data points an in- 
clusion of smaller scales will enhance constraints from {'R'Re) 
stronger than those from (M%„). 

The iiB-statistics gives tighter constraints on cosmology than 
the optimized ring statistics, which can be explained by its 



Table 1. Values of q resulting from the likelihood analyses of 
the 5 data vectors. 



Data vector q 

(•KRe) (t? = %^ = 1') 1518 

(KRe) (r, = iWY, 'P™ = 6') 177.3 

(KKe) (tl = 0-1) 207.9 

(M%) (6 mm = 6') 169.8 

E CP,™ = 1') 122.5 



larger S/N ratio. However, we do not use the EB-statistics to 
analyze the CFHTLS data in the next section for the reason 
that the /sB-statistics' data points are strongly correlated (see 
Fig. |3). In order to identify B-modes as a function of angular 
scale accurately, the lower correlation of the ring statistics is 
more useful. 

7. Ring statistics with the CFHTLS 

In Sect. 15. II we have shown that the ring statistics' data points 
are significantly less correlated compared to data points of the 
aperture mass dispersion. Therefore, despite its lower S/N ra- 
tio, the ring statistics provides an ideal tool to analyse B-mode 
contaminations depending on the angular scale. In this section 
we use the 2PCFs of the FSH08 analysis and therefrom calcu- 
late the ring statistics for a scale-dependent 77 = ^min/^ and 
for 77 = 0.1. We performed a similar analysis for other cases of 
fixed 77, which resulted in a significantly weaker signal. 
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Fig. 6. The ring statistics signal measured from the CFHTLS for the case of rj — # m in/ v I' {upper row). The red data points (circles) 
correspond to the E-mode signal, the black data points (triangles) to the B-mode signal. The three panels correspond to small 
{left), intermediate {middle), and large {right) scales. The lower row shows a similar analysis but for 77 = 0.1. 



The CFHTLS 2PCF was measured in 72000 bins over an an- 
gular range of 0:05 < # < 466'; we calculate {KKe) (Eq. 
17) and {KKb) (Eq. [8} in 60 logarithmic bins over a range 
0f5 < ¥ < 460'0. The error for the ;'-th E/B-mode data point is 
calculated as ^/Cr e/b (*P,', *?,), where Cr e/b ( ( P,', SP,-) is calculated 
from a Gaussian 2PCF covariance. This Gaussian covariance 
was calculated from a theoretical model using the same cos- 
mology and survey parameters as in the FS H08 analysis. We 



do not employ the non-Gaussian correction of lSemboloni et al 
(12007b as this corrects the C++ -term in the 2PCF covariance, 

but not the C - and C+_-terms. Here, we use the full 2PCF 

covariance in the analysis. Similar to FSH08 we do not con- 
sider systematic errors in our analysis which might lead to an 
underestimation of the error bars. 

The results of our analysis are illustrated in Fig. [6] The three 
panels in the upper row show the ring statistics' E- and B- 
modes on (from left to right) small, intermediate and large 
scales of VP for the case of rj — J? m i n / V f. The three panels in 
the lower row show the same analysis but for 77 = 0.1. The 
circled (red) data points correspond to the E-mode signal, the 
triangled (black) data points correspond to the B-mode signal. 
We measure a robust E-mode shear signal, however we also 
find a significant B-mode contribution on small (around 2'), in- 
termediate (16'-22'), and large scales (right panel). On small 



scales E-and B-mode are of similar order. It should be stressed 
that such an analysis of small-scale contaminations is not feasi- 
ble with the aperture mass dispersion, which, to avoid the E/B- 
mode mixing on small scales, involves a theoretical (therefore 
B-mode free) 2PCF in its calculation. This theoretical data ex- 
tension, combined with the fact that the aperture mass disper- 
sion data points are stronger correlated (Sect. [5]) can hide pos- 
sible small-scale contaminations in the data. 
The B-mode contamination on large scales is also observed in 
the FSH08 analysis. In addition, we find a small B-mode on 
intermediate scales (between 16' and 22'), otherwise these in- 
termediate scales are mostly free of B-modes and give a ro- 
bust E-mode signal. The small correlation of the individual data 
points leads to the oscillations in the amplitude of the shear sig- 
nal. A similar analysis with the aperture mass dispersion shows 
a much smoother behavior. 



8. Conclusions 

Decomposing the shear field into E- and B-modes is an impor- 
tant check for systematics in a cosmic shear analysis. The most 
commonly used methods for E- and B-mode decomposition, 
namely the aperture mass dispersion and the E/B-mode shear 
correlation function, require the 2PCF to be known down to 
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arbitrary small or up to arbitrary large angular separations. 
In practice, the 2PCF is only measured over a finite interval 
[)? m i n ; #max]- As a result the aforementioned methods do not 
separate E- and B-modes properly, e.g. the aperture mass 
dispersion suffers from E/B-mode mixing on small angular 
scales (see KSE06 for further details). 

In contrast, the ring statistics (invented in SK07) separates 
E- and B-modes properly using 2PCFs measured on a finite 
angular scale. As outlined in SK07 the filter functions of the 
ring statistics, i.e. Z ± , are in general complicated to calculate; 
the authors restrict the free parameters this filter function to 
one, namely 77. This parameter is held fixed, independent of the 
angular scale *F at which the ring statistics is evaluated. In this 
paper, we improve on the condition of a fixed 77 by choosing a 
scale-dependent 77 = # m j n / v F which significantly improves on 
the ring statistics' S/N ratio. 

Furthermore, we present a formula to calculate the ring 
statistics' covariance from a 2PCF covariance. This formula 
is applied to a 2PCF covariance obtained from ray-tracing 
simulations. We therefrom calculate the correlation matrices 
of ring statistics and aperture mass dispersion and find that 
the data points of the first are significantly less correlated than 
the data points of the second. We employ these covariances 
to compare the information content of the two second-order 
statistics and find that the ring statistics' data points place 
tighter constraints on cosmological parameters than data points 
of the aperture mass dispersion. The reason for this is that we 
can include smaller scales in the ring statistics's data vector 
which is not possible for (M% p ) due to the aforementioned 
E/B-mode mixing. In addition, we consider a polynomial 
filter function which decomposes E- and B-modes on a finite 
interval and therefrom calculate an additional second-order 
measure, the EB-statistics. We compare the correlation of 
data points and the information content of this EB-statistics 
to the ring statistics and find that it shows a significantly 
larger correlation of the data points, but a higher information 
content. This can be explained by the high S/N ratio of the 
.EB-statistics. 

We apply the ring statistics with 77 = J? m i n /T' and 77 = 0. 1 
to CFHTLS data, more precisely we calculate both from the 
2PCF used in the latest CFHTLS analysis (FSH08). We mea- 
sure a clear shear signal for 77 = J^min/*? which decreases when 
performing the same analysis for 77 = 0.1. The fact, that data 
points of the ring statistics have small correlations enables us 
to determine the contaminated scales very accurately. We find 
B-modes on large scales which is comparable to the findings of 
FSH08. In addition, we detect B-modes on intermediate (16' 
- 22') scales and a scattered B-mode contribution on scales 
below 3'. In the latter case the shear signal is of the same order 
as the B-mode contribution. 

A similar analysis with the aperture mass dispersion is only 
possible when including a 2PCF from a theoretical model 
in order to avoid the E/B-mode mixing on small angular 
scales. These added theoretical data can conceal remaining 
systematics (B-modes) which can be identified properly using 
the ring statistics. This property is most likely the most useful 
feature of the ring statistics. It can be used to detect remaining 
systemics very accurately in future surveys. 



The noise-level of the ring statistics on small scales can be 
reduced by increasing the number of galaxy pairs within the 
contributing 2PCF-bins. The number of galaxy pairs inside a 
2PCF-bin increases quadratically with « ga i, therefore it would 
be interesting to test the ring statistics on a data set like e.g. 
the COSMOS survey. Similarly, an increased survey volume 
will significantly enhance the constraints, for the reason 
that the cosmic variance scales with I /A. For example, the 
CFHTLS data we consider here covers an area of 34.2 deg 2 
with 77 ga i = 13.3. Testing the ring statistics on the full CFHTLS 
sample (172 deg 2 ) would be an interesting project in the future. 



Appendix A: T+ -functions 

In order to define the r ± -functions used for the calculation of 
the .Efi-statistics we remap 1? 6 [i? m j n ; i? max ] to the x e [-1; 1] 
and define 

2 ft - ft ■ -ft 

z. u f^min '-'max 



B = 



ft - ft ■ 

"max [/ min 

ft -ft ■ 
u max u min 

$max "I" ^min 



(A.l) 
(A.2) 



We choose our filter function T + (x) to be the lowest order poly- 
nomial which fulfills the two integral constraints of Eq. ( TT3l 
and the normalization J dx T + (x) T + (x) = 2. The function 
reads 

T + (x) = (3B 2 -5-6Bx + 3(5- B 2 )x 2 ) , (A.3) 



where 



! (25 + 5 B 2 + 6 Br) 



(A.4) 



Given the analytic form of T+ the corresponding T is uniquely 
determined through Eq. 0141) . 
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